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Abstract 

Using lattice effective field theory, we study the ground state binding energy of N distinct 
particles in two dimensions with equal mass interacting weakly via an attractive (iV)-symmetric 
short range potential. We find that in the limit of zero range and large N, the ratio of binding 
energies Bn/Bn_i approaches the value 8.3(6). 



1 



I. INTRODUCTION 



We consider the ground state of N distinct particles in two dimensions with equal mass 
interacting weakly via an attractive SU (N)-symmetTic short range potential. Since the 
ground state is completely symmetric this is equivalent to the problem of N weakly-bound 
identical bosons. The self-bound two-dimensional system with a realistic van der Waals 
potential is relevant to the properties of adatoms on surfaces. In this work, however, we 
address the question of what happens as the range of the interaction goes to zero, 

V(x ir -- ,x N )^C ^(xi-xj). (1) 

l<i<j<N 

Let -Bat be the ground state binding energy of the iV-particle system in the zero range 
limit. The first calculation of B 3 /B 2 was given in [l | . The precision of this calculation was 
improved by and most recently a precise value of B3/B2 = 16. 522688(1) was given in 
0| . There have also been studies of the four- and five-particle systems j^O,!^. But range 
corrections for these studies appear to be very large, and the first precise determination of 
B^/B 2 in the zero range limit was only recently given in jjj, yielding a value of B±/B 2 = 
197.3(1). 

n 

The behavior of Bn in the large- N limit was also recently discussed in |3J. They showed 
that due to the weakening of the attractive coupling at short distance scales, the large- N 
droplet system could be treated classically. This yielded a prediction for the ratio of the 
binding energies in the large- N limit, 
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lim — — ~ 8.567. (2) 

A^oc B N _i 



the iV-particle system for N < 7 was investigated using diffusion Monte Carlo with 
both a Lennard- Jones potential and a more realistic helium-helium potential. However 
the results showed that range corrections were too large to allow for a determination of 
B N /B N _x for large N. 

Although there is no known system of atomic or molecular clusters that displays the 
physics of the zero range limit for large N, the topic is interesting for several reasons. 
With recent advances in laser trapping techniques it is now possible to produce many-body 
quantum systems on a two-dimensional optical lattice. Much of the attention has been 
devoted to the Bose-Hubbard model with repulsive on-site interactions 0, [ljj], but the 



weakly attractive iV-boson system can also be studied. In that case computational lattice 
studies such as this would be of immediate relevance. Our system also raises interesting 
questions about the convergence of effective field theory and the large- N limit. Results of 
previous numerical studies suggest that it is surprisingly difficult to reach the zero range and 
large- N limits at the same time. We explore why this is the case and what can be done to 
overcome some of the difficulties. Similar issues arise in systems of higher-spin fermions in 
optical traps and lattices. In these systems the competition between short range interactions 
and large- N effects can determine properties of the ground state, two-particle pairing versus 
multi-particle clustering [ill Q| • 

In this paper we study the iV-particle system using lattice effective field theory. The 
organization of our paper is as follows. We first discuss the renormalization of the interaction 
coefficient in the two-particle system. We discuss renormalization in the continuum with 
a sharp momentum cutoff and then on the lattice. After that we address two features of 
the large- N limit. The first is a rescaling technique that cancels some of the nonzero range 
corrections from the ratio / B^-i- The other is an overlapping interaction problem 
that occurs when many particles lie within a region the size of the range of the potential. 
We show that this problem can produce large systematic errors that grow with N. The 
strength of the overlapping interaction must be reduced if we wish to probe zero range 
physics accurately for large N. We demonstrate one way of doing this which exploits an 
unusual feature of the discrete Hubbard-Stratonovich transformation Q]. Using lowest- 
order lattice effective field theory, we compute B^ / Bn_i for N < 10. Extrapolating to the 
limit N — > oo, we find the result 

lim -^L = 8 .3(6). (3) 

II. TWO-PARTICLE SYSTEM AND RENORMALIZED COUPLING 

We begin by reviewing the two-particle system in the continuum formalism with a sharp 
cutoff, A, on the magnitude of the momentum. For a zero range potential, 

V(x ir -- ,x N )=C ^(xi-xj), (4) 

l<i<j<N 

the diagrams which contribute to two-particle scattering are shown in Fig. ^ We let m be 
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FIG. 1: Diagrams contributing to two-particle scattering for a zero-range potential. 



the particle mass. In order that the bound state pole in the rest frame occurs at energy 
E = —B 2 , we get the constraint 
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We can rewrite this as 



mB 2 A" 



1 + mB 2 A- 2 
and the bound state energy is given by 
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We now consider the same calculation on the lattice. Let a be the spatial lattice spacing 
and a t be the temporal lattice spacing. We start with the Hamiltonian lattice formulation 
where a t — 0. The standard lattice Hamiltonian with nearest neighbor hopping has the 
form 

n l<i<N l=x,y 

+ Ca ' 2 H E &!(*om*o&5(*)W (8) 

n l<«<j<A r 

Here, bi(n) is an annihilation operator for a particle with flavor i at the spatial lattice site 
n. The condition on C in the Hamiltonian lattice formalism is 



— y — 

2r,2 R„ i ' 



— = lim 

C l-oo a 2 L 2 , ^ 5 2 + 2fi_ ' 

k integer 



(9) 



where fL is the lattice kinetic energy and L is the length of the periodic lattice cube in 



lattice units. For the standard lattice action 



k ma 2 ^ L L J 



(10) 



s=x,y 



For later reference we define uo as the momentum independent term inside the summation 
in (jlOj) . So for the standard lattice action u — 1. We define the lattice cutoff momentum 
A = 7ra _1 . Then in the limit A — > oo, 



mi^A 2 = B exp 
where for the standard action B ~ 3.24. 
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(m£ 2 A- 2 ) 2 , (11) 



In order to test the cutoff dependence of our lattice results, we also consider actions 
with 0(a 2 )-improved and 0( a 4 ) -improved kinetic energies. 0(a 2 ) and 0(a 4 ) corrections 
to the interaction are not included since these would entail a significant number of new 
interactions. Because we are not performing a full 0(a 2 ) or 0(a 4 ) improvement, we do 
not expect the improved kinetic energy actions to give qualitatively better results than the 
standard action. However a comparison of the different actions provides an additional check 
that the results reproduce continuum limit behavior rather than lattice-dependent artifacts. 
For the 0(a 2 )-improved action the lattice kinetic energy is 

«'4E [f-loos^ + icos^]. (12) 



k 



=x,y 

In this case uj = | and B ~ 1.79, where B is defined in the asymptotic expression (J11J) . For 
the 0(a 4 )-improved action 

^: = ^E [§-§ cos ^ + i cos ^-MCOS^], (13) 

s=x,y 

ou = g|, and B ~ 1.54. As we increase the order of improvement, fL more closely approxi- 
mates the continuum kinetic energy and B approaches the continuum sharp cutoff value of 
1. 

At nonzero temporal lattice spacing the same diagrams in Fig. ^ contribute to two- 
particle scattering. A derivation of the Feynman rules at nonzero temporal lattice spacing 
for the analogous three-dimensional system can be found in jl4| ]. as well as a derivation of 
the bound state pole condition. In two dimensions the strength of the interaction is given 
by the transfer matrix element 

e-^- 2 -l)(l-u^)\ (14) 
/ \ ma z J 

while the free lattice propagator has the form 



L t is the total number of temporal lattice units and ko is an integer from to L t — 1. 

As we take L t — > oo, the energy in physical units becomes a continuous variable. Re- 
quiring that the bound state pole in the rest frame occurs at energy E = —B 2 , we get the 
constraint 
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At nonzero temporal lattice spacing we therefore have 
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where 
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In this case B is a function of atm _1 a -2 and is different for the standard and improved 
lattice actions. For given values of B 2 ,a, and at, we determine C in the infinite volume 
limit L — ► oo. For mB 2 hr 2 not too small, roughly 10~ 6 or larger, we use the exact 
expression (fTB^) for sufficiently large values for L. For smaller values of mB 2 A~ 2 it is more 
convenient to use the asymptotic expression ()17j) . But once the interaction coefficient C 
is determined and we proceed to the iV-body system, such extremely large lattice volumes 
are unnecessary. It suffices to consider lattice systems larger than the characteristic size 
of the iV-body droplet. For large iV this is many orders of magnitude smaller than the 
characteristic size of the two-body droplet. 



III. RATIOS IN THE LARGE-iV LIMIT 



It has been suggested that the large-A ground state wavefunction can be described as a 
dassica! d^ibution & If « N is the c_ ist ic S i Z e of the droplet, the d^ibution , 
proportional to ip(r/R^) for some function if) independent of N, and the binding energy B^ 
is proportional m~ 1 R^ f 2 . In order to determine ip, one integrates out high energy modes 
to determine the effective coupling at energy B N . If this picture of the large- A droplet is 
correct, then errors due to the finite cutoff momentum A appear only in the combination 
mi^jvA^ 2 . Therefore if we measure binding energies while keeping mi^ArA -2 fixed, much of 
the error cancels in the ratio B^ / B^_i. In essence we are using large- A similarity under 
rescaling to eliminate cutoff errors. If the classical droplet picture is incorrect, then this 
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technique will probably not reduce errors. The issue will be settled when we analyze results 
of the Monte Carlo simulations. 

Let -Bjv(A) be the measured binding energy of the iV-particle ground state at cutoff 
momentum A. Conceptually it is simplest to regard m and B 2 as fixed quantities while we 
vary A. In the continuum limit 

lim B N (A) = B N . (19) 

A^oo 

Let z > be a parameter that measures proximity to the continuum limit, 

z = mB N (A) ■ A' 2 . (20) 
For a given z, we define the cutoff momentum A(z, N) implicitly so that 

mB N (A(z, AO) • (A(z, N))' 2 = z. (21) 

We define f(z) as 

f(z) = lim i In [B N (A(z, N))/B N ] . (22) 

N^oo iV 

f(z) measures the exponential growth of finite cutoff errors with increasing N. We have 



lim < In 



B N {A{z,N)) 



and so 



lim 



B N 
B N 
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e~ f{z) lim 



Bn-i 
B N (A(z,N)) 



(23) 



(24) 



n-^o B N -i n-kx B N -i(A(z,N — 1))' 

Therefore so long as \f(z)\ 1, the large- N ratio of binding energies can be measured 

reliably. Other cutoff errors which do not grow linearly with iV will cancel in the ratio 

B n ^(A( Zj N-1)) { ] 

as we take iV — > oo. 

In our Monte Carlo lattice simulations it is more convenient to regard m and A as fixed 
quantities while varying B 2 . We define B 2 (z,N) implicitly by 

mB N (B 2 (z,N)) ■ A' 2 = z. (26) 

We are changing the overall physical scale when we change B 2 , and so we work with ratios 
BnI ' B 2 . The analog of the result is 

B N __ fM ^ B N (B 2 (z,N))/B 2 (z,N) 



lim 

TV^oo Bn-1 



e" /(2) lim 



iv^oo B N ^(B 2 (z, N - 1))/B 2 (z, N-V 



(27) 
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IV. OVERLAPPING RANGE AND IMPLICIT iV-BODY INTERACTION 



Large range corrections can occur when many particles lie within a region the size of 
the range of the potential, A -1 . The problem is most severe when all iV particles lie in 
this localized region, and the potential energy is amplified by a factor of N(N — l)/2. For 
a continuum potential with a repulsive core, the result is a deep hole at the center of the 
multiparticle wavefunction and a tendency towards underbinding or unbinding. At lowest 
order in lattice effective field theory the effect goes in the opposite direction. A spike forms 
at the center of the wavefunction when all particles lie on the same lattice site, and the 
binding energy is too large. 

Consider the state with A" particles at the same lattice site in the Hamiltonian lattice 
formalism, 

\U N ) = b\(n)bl(n)--.bUn)\0). (28) 
The expectation value of the potential energy for this localized state is 

(n »|y|n»> = ™^. (29) 

This potential energy can be regarded as an implicit A^-body contact interaction produced 
by overlapping two-body interactions. In the continuum limit we know that the importance 
of this A^-body contact interaction is suppressed by many powers of the small parameter 
z = mB N A~ 2 . However the situation at finite A can be quite different from the continuum 
limit if the potential energy per particle for the localized state (n^) is as large as the cutoff 
energy scale, 



C(N -I) 



2a 2 

To lowest order in mB 2 A~ 2 , the renormalized coupling is 

An 

C 



(30) 



mln(mB 2 A 2 ) 

An 



mm (itiBnA^ 2 ) — mm (Bn/B 2 ) 



(31) 



For large A" 



where 



An 

(j ~ ; (32) 
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(3 = Mm p-. (33) 
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Then 



C(N-l) 



7T 



2o? ma 2 
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ln/3 



In the continuum limit the problem goes away since 

1 



0. 



(34) 



(35) 



ln/3 — jj: In z 

However the convergence is slow and requires z <C e~ N . For actual lattice simulations it is 
therefore necessary to limit the size of the implicit iV-body contact interaction. 



V. DISCRETE HUBBARD-STRATONOVICH TRANSFORMATION 



There are several ways to deal with the large implicit iV-body contact interaction. On 
the lattice there is one method which is particularly convenient. This is to write the two- 
body interaction using a discrete Hubbard- Stratonovich transformation 13j. The discrete 
Hubbard-Stratonovich reproduces the two-body contact interaction exactly. Typically it 
is used for systems with spin- 1/2 fermions where Pauli exclusion implies that there are no 
iV-body contact interactions beyond N = 2. It seems therefore that the properties of 
the transformation for N > 3 has not been discussed in the literature. In the following 
we show that when a discrete Hubbard-Stratonovich transformation is used, the temporal 
lattice spacing regulates the strength of the implicit iV-body contact interaction. 

For simplicity we show only the interaction part of the Hamiltonian. The exponential of 
the two-body interaction at site n over a Euclidean time step a t is 



e -atH int _ exp 



—a t Ca 2 b\{n)bi{n)b^{n)bj(n) 

l<i<j<N 



(36) 



The discrete Hubbard-Stratonovich transformation amounts to making the replacement 



-<hH- w 



\ J2 ex P " (\a t Ca- 2 + \s(n) \ I £ &J(*0&i(*0 ~ 1 ) 

s(n)=±l L ^ ' \l<i<N J 



where 



cosh A ox]) ( —-a t Ca 2 



A > 0. 



To see that this has all the desired properties, let us define 



MK)= 1 - £ 



exp 



s(n)=±l 



-a t Ca- 2 + Xs(n) ) (K - 1) 



(37) 



(38) 



(39) 



for nonnegative integer K. We note that A(0) = A(l) = 1, and A(2) = exp (—atCa 2 ). 
These are precisely the expectation values of e~ atHint for K = 0,1,2 distinct particles at 
lattice site n. When K > 3 but \{K - 1) < 1, we find 



A(K) ~ exp 



-a t Ca 



_ 2 K(K-1) 



(40) 



This is also the expectation value of e atHint for K distinct particles at lattice site n. However 
when K > 3 and \{K - 1) > 1, 



A(K) ~ - exp 



-a t Car 2 + A ){K - 1 



(41) 



The total potential energy of the i^-particle localized state, jn^, no longer increases 
quadratically with K. The temporal lattice spacing a t acts as an auxiliary ultraviolet 
regulator that limits the size of the implicit ii'-body contact interaction. When K < 2 or 
the implicit i^-body contact interaction is smaller than a^ 1 , we have the unaltered result, 

CK(K - 1) 



{li K \V\ll K ) 



2a? 



(42) 



When K > 2 and the implicit i^-body contact interaction exceeds a t , then the regulator 
takes effect and we have 



(T1 K \V\T1 K ) 



1 



-atCa 



A (K- l)+ln2 



(43) 



VI. ALGORITHM 



The standard lattice action we use for our simulations is 



E 



n,nt,i 



c*(n,n t )ci(n,n t + 1) - e 



^t^-Xs{n,n t ) ( 1 2,/, 



2 I c^ntjCi^nt) 



at 



2ma 2 



\c*(n, n t )ci(n + 1, n t ) + c*(n, n t )ci(n — I, n t )j — ^ As(n, n t ), (44) 



n,nt,l,i 



n,n t 



where rif is the temporal lattice coordinate, Cj is the path integration field for the particle 
of type and s is the discrete Hubbard-Stratonovich field which takes values ±1. We have 
used the ,a« C e conveys d eve lo ped in QQfora th.ee-di.nensional lattice. The ehoiee 
of Bose/Fermi statistics for q is irrelevant since we consider systems with no more than one 
particle of each type. 
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In order to compute the ground state binding energy Bn we consider the correlation 
function 

Z N (t) = (*° N \ e~ Ht \K) , (45) 
where the initial/final state is the state with all N particles at zero momentum, 

|*Sr> = &I(0)S(0) • • • 6^(0)10). (46) 

\*&%) is also the ground state of the non-interacting system. We refer to t as Euclidean time 
and define 

E N (t) = -^[\nZ N (t)}. (47) 

Then as t — > +oo, E^(t) converges to —-Bat, the ground state energy of the interacting N- 
particle system. The only assumption is that the ground state has a nonvanishing overlap 
with the ground state of the non-interacting system. 

The conversion of t 
has been discussed in 

we keep the discrete Hubbard-Stratonovich field and calculate the sum over configurations, 



re lattice action to a transfer matrix formalism at fixed particle number 
16 . We use the same transfer matrix derived there, except in this case 



Z N (t) oc^e-^»* As( "> nt) (^|r[e-^] \V° N ), (48) 

s 

H(s) consists of only single-body operators interacting with the background Hubbard- 
Stratonovich field. We can write the full iV-particle matrix element as the N th power 
of the single-particle matrix element, 

(*° N \T[e- H W}\*° N ) oc [M(s,t)] N , (49) 



M(s,t) 



k = 



T [e- H ^ 1 ] 



k = 



(50) 



where k = Oy is a single-particle state with zero momentum. Our time-ordered exponen- 
tial notation, T [e - ^^*] , is shorthand for the time-ordered product of single-body transfer 
matrices at each time step, 



T [e- H ^] = M (Lf _D ■ . • • • M (nt) • . . . • M (1) • M (0) , 



(51) 



where L t is the total number of lattice time steps and t = L t a t . If the particle stays at the 
same spatial lattice site from time step n t to n t + 1, then the corresponding matrix element 
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of M( nt) is 



atCa . , / 2a+ 
- — As(ra,nt) I i - < 



ma 2 



e -^^-Xs { n,n t) I 1 _L j (52) 



If the particle hops to a neighboring lattice site from time step to n t + 1 then the corre- 
sponding matrix element of Mt nt ) is j^ir- All other elements of M(„ t ) are zero. 
We sample configurations according to the weight 

expi J2 Xs (^ n t)+ Nl °s[M(s,t cnd )} I, (53) 

^ ra,n t J 

where t en d is the largest Euclidean time at which we wish to measure Z^(t). We use a 
simple heat bath/Metropolis update procedure. For each configuration the observable that 
we compute is 

Q(.,t)= |M(M>r jV . (54) 

for t < t en( j. This is the same general technique that was used in |l7[. By taking the 
ensemble average of 0(s, t) we are able to calculate 

Z N (t) 



ZN(t en d) 



(55) 



VII. RESULTS 



For each simulation we have computed roughly 2 x 10 5 successful heat bath/Metropolis 
updates for each lattice site, split across four processors running completely independent 
trajectories. Averages and errors were calculated by comparing the results of each processor. 
The codes were based on existing codes used for light nuclei in three-dimensions and we have 
kept some of the same input parameters relevant for the light nuclei system. We use a mass 
of m — 939 MeV and keep the spatial lattice spacing fixed at a = (40 MeV) -1 . This 
corresponds with A = ira^ 1 ~ 126 MeV and cutoff energy A 2 /m = 16.8 MeV. Clearly 
these input parameters in raw form are not appropriate for atomic clusters. Therefore we 
translate of all of the parameters in terms of dimensionless numbers which can then be easily 
applied to any two-dimensional droplet system. 

We have already defined the dimensionless ratio z, 

B^ma 1 Tx^' 1 . (56) 



A 2 /m 
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z measures the ratio of Bn to the cutoff energy. In most cases it is clear which iV we are 
referring to and so we use the simple notation z. When there is some possibility of confusion 
we include the N subscript, z^. 
We also define e, 



e = ™ l = a^ma 2 ^' 1 . (57) 



A 2 /m 

A small value for e indicates that there is a significant amount of high frequency regulariza- 
tion provided by the nonzero temporal lattice spacing a t . A large value for e means that 
we are close to the Hamiltonian limit, a t — > 0. There is little or no regularization of high 
frequency modes and most of the regularization is provided by the momentum cutoff A. 

We adjust the two-particle binding energy B 2 in order to study the finite cutoff depen- 
dence. Since we keep A fixed, our value B 2 will decrease as go to larger values of N. For 
convenience we use the shorthand 

b N = B N /B 2 (58) 

for the dimensionless ratio of the binding energies. For each data point we increase the 
spatial length and temporal extent of the lattice until the finite volume/time errors are clearly 
smaller than the statistical errors. The largest lattice system we simulate is 9 x 9 x 260. 

We have computed for A/" < 10 for a wide range of values for B 2 using the 0(a 4 )- 
improved action and a t = (20 MeV) -1 , which corresponds with e = 3.7. The results are 
shown as a plot of ln(b^) versus z in Fig. El We see that there is considerable dependence 
on z. The dependence appears to be roughly linear in z for 0.1 < z < 0.3, and we have 
drawn interpolating lines. We note that since ln(6jv) and ln(fojv-i) have approximately the 
same slope, most of the z dependence cancels in the combination ln(fojv) — hi(&7v-i)- This 
suggests that / (z) as defined in ()22|) is small. Much of the systematic cutoff errors can be 
cancelled in the ratio ^/^iv-iby keeping z the same for b^ and 6jv-i- From Fig. |2]we see 
that b^/b^-i is about 10 for 5 < A^ < 10. Therefore scaling B 2 proportional to 10 _Ar as 
we probe the A^-body droplet should keep z approximately the same for these values of N. 

Next we calculated 6jv/fcjv-i f° r A^ < 10 using three different actions. We compared the 
standard action, the 0(a 2 )-improved action, and the 0( a 4 ) -improved action, using a t = (20 
MeV)" 1 and B 2 = 2 x lO 2 "^ MeV. This corresponds with e = 3.7 and z 2 = 1.2 x lO 1 ^. 
The results are shown in Fig. El We see about a 10% variation among the three different 
actions, with the 0(a 2 )- and 0(a 4 )-improved actions agreeing slightly better with each other 
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FIG. 2: ln(b_y) versus z. We use the 0(a 4 )-improved action and e = 3.7. 



than with the standard action. 

In Fig. HJwe plot 6_v/6_v-i using the 0{ a 4 ) -improved action, a t = (20 MeV) -1 , and three 
different sets of values for B 2 : B 2 = 3 x lO 2 "^ MeV, 2 x 10 2 ^ MeV, and 1 x lO 2 ^ 
MeV. This corresponds with e = 3.7 and z 2 = 1.8 x lO 1 -^, 1.2 x lO 1 -^, and 0.6 x lO 1 ^ 
respectively. The discrepancies for the different values of B 2 are at the 30% level for small 
N, but as expected the errors decrease with increasing N. 

We also studied the dependence of &_v/&jv_i on the temporal lattice spacing a t . We 
set B 2 = 2 x lO 2 "^ MeV and used the 0(a 4 )-improved action with a t = (16 MeV)- 1 , (20 
MeV)" 1 , (30 MeV)" 1 , and (40 MeV)- 1 . This corresponds with z 2 = 1.2 x lO 1 -^ and 
e = 3.0, 3.7, 5.6, and 7.5 respectively. Since e is rather large, a. has only a small effect on 
the ultraviolet regularization of the two-body interaction. Instead the importance of at is 
as an auxiliary regulator on the implicit iV-body contact interaction. The results are shown 
in Fig. E| The results appear to differ at about the 10 — 15% level. 

In Fig. |H1 we combine all of the data shown in Figs. 01 EJ and 03 For comparison we 
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FIG. 3: Comparsion of /bjv-i for the standard, 0(a 2 )-improved, and 0(o 4 )-improved actions. 
We use e = 3.7 and z 2 = 1.2 x lO 1 ^. 



22 
20 
18 
16 - 
14 - 
12 



> 10 



8 - 
6 - 
4 



i 1 1 1 r 

z 2 = 1.8 x 10^ i — h- 
4 = 1.2x10 "J H-x- 
z 2 = 0.6x10 i yv :— *- 



t r 



3E 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

MN 



FIG. 4: Comparsion of b^/hx-i for different values of 22- We use the 0(a 4 )-improved action and 
e = 3.7. 
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FIG. 5: Comparsion of b^/b^-i for different values of e. We use the 0(a 4 )-improved action and 
z 2 = 1.2 x lO 1 ^. 



include the known results for N = and N —> oo [31. We draw two best fit 

curves with up to quadratic dependence on 1/N. The known results were not included in 
this fit. The best fit curve using 1/N and 1/N 2 gives a value 



hm ~ 



7.7, 



while the best fit curve using only 1/N 2 gives a value 

lim ^ N ~ 8.8. 

If we take these two results as approximate lower and upper bounds then we find 



hm ~ 

N-><x> On-1 



5.3(6). 



(59) 



(60) 



(61) 



VIII. CONCLUSIONS 



We have studied the two-dimensional iV-particle system with short range attraction using 
lowest-order lattice effective field theory. We discussed two aspects of the large- N limit. The 
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FIG. 6: Cumulative data for 6^/^-1- We draw two best fit curves with up to quadratic depen- 
dence on 1/N and show the known results for N = 3, 4, and oo. 



first is a technique that uses large- N similarity under rescaling to cancel some of the nonzero 
range corrections from the ratio B^/B^-i- The other is the problem of a large implicit 
iV-body contact interaction when many particles lie within a region the size of the range of 
the potential. We regulated this implicit iV-body contact interaction on the lattice using 
a discrete Hubbard-Stratonovich transformation. Using a heat bath/Metropolis algorithm 
we computed Bn/Bn-i for N < 10. Extrapolating to the large- N limit we found 

lim — — = 8.3(6). (62) 

N^oo B N -! 

This appears to be in agreement with the value 8.567 found by Q|. 

While we have measured the large- N limit of B^ / B^-i to within 10%, we relied on large- 
N similarity under rescaling to keep the finite cutoff errors in check. The z dependence 
in Fig. |21 suggests that one needs to go beyond leading order to accurately describe all of 
the physics at large N. This competition between effective field theory expansions and 
the large- N limit presents an interesting theoretical challenge. Since there are no known 
physical systems where we can experimentally measure the universal zero range behavior, 
the coefficients of the higher- dimensional operators must be set by numerical calculations. 
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One technique perhaps is to use numerical renormalization group matching to relate the 
coefficients of higher-dimensional operators for different values of m£> 2 A~ 2 . However more 
study would be needed to see if this is a viable technique. 
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